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Abstract 

We obtain free of resonances regions for the elasticity system in the exterior of a strictly 
convex body in R 3 with dissipative boundary conditions under some natural assumptions 
on the behaviour of the geodesies on the boundary. To do so, we use the properties of the 
parametrix of the Neumann operator constructed in [T2]. As a consequence, we obtain time 
decay estimates for the local energy of the solutions of the corresponding mixed boundary 
value problems. 

1 Introduction and statement of results 

Let O C R 3 be a strictly convex compact set with smooth boundary V = dO and denote 
by Q = R 3 \ O the exterior domain. Denote by A e the elasticity operator, which is a 3 x 3 
matrix-valued differential operator defined by 

A e u = fj> Au + (A + /U )V(V • u), 

u — (ui, 1*2,1*3). Here Ao, fJ-o are the Lame constants supposed to satisfy 

Ho > 0, 3A + 2/i > 0. (1.1) 

The Neumann boundary conditions for A e are of the form 

3 

{Bu)i\ T :=£>ii(«K-|r = 0, i = 1,2,3, (1.2) 
3=1 

where 

, . . _ . / dui duj\ 

is the stress tensor, v is the outer unit normal to T. The purpose of the present paper is to study 
the time decay properties of the elasticity system in with dissipative boundary conditions. 
More precisely, we are going to study the following mixed boundary value problem 

{pi - A e )u = in (0, +00) x ft, 

Bu-iAu = on (0, +00) x T, (1.3) 
u(0) = /1, d t u(0) = f 2 , 
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where A is a classical zero order 3x3 matrix- valued pseudo-differential operator on T, indepen- 
dent of t and satisfying the properties A = A* , A > 0. Moreover, we suppose that there exist a 
non-empty compact set Tq C T and a constant C > so that we have 



(A/,/) L2(r) >C||/|| 2 L2(ro) . 



(1.4) 



The large time behaviour of the solutions to (1.3) with A = is well understood. Kawashita 
[6] showed that there is no uniform local energy, while Stefanov-Vodev [12], [13] proved the 
existence of infinitely many resonances converging polynomially fast to the real axis. The reason 
for this is the existence of surface waves (called Rayleigh waves), that is, a propagation of 
singularities of the solutions along the geodesies on T with a speed cr > strictly less than 
the two other speeds in S7. Therefore, a strictly convex obstacle is trapping for the Neumann 
problem of the elasticity wave equation. Note that for the Dirichlet problem it is non-trapping, 
and in particular we have an exponential decay of the local energy similarly to the classical wave 
equation (see [17J). Comming back to the equation (1.3) with non-trivial A, note that we still 
have a propagation of singularities of the solutions along the geodesies on T with a speed cr > 0. 
Therefore, in order to be able to get a better decay of the local energy we need to suppose that 
all geodesies meet the part on T where the dissipative term is non-trivial. More precisely, we 
suppose that there exist a non-empty open domain r' C Tq and a constant T > so that 

for every geodesies 7 with 7(0) G I\ there exists < t < T such that j(t) G r' . (1.5) 

The outgoing resolvent, R(X), corresponding to the problem (1.3) is defined via the equation 



Recall that "A-outgoing" means that there exist a 1 and a compactly supported function g 
so that 

R Wf\\x\>a = RoWg\\x\> a , 
where i?o(A) is the outgoing free resolvent, i.e. 

R (X) = (A e + A 2 )" 1 G C{L 2 ) for ImA < 0. 

Let x £ Co^R- 3 )) x = 1 on C. in t ne same way as in the case A = we have that the cutoff 
resolvent 



extends meromorphically to the whole complex plane C with poles in Im A > called resonances. 
One of our goals in the present paper is to study the distributions of the resonances near the 
real axis under the above assumptions. Our first result is the following 

Theorem 1.1 Under the assumptions (1-1), (1-4) an d (1-5), Rxity extends analytically to 
{|ImA| < Ci|A| _1 , |ReA| > C2 > 0} and satisfies there the estimate 




(1.6) 



R X (X) := X R{X)x 



(1.7) 
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Moreover, under the assumption (1.1) only, there exists a constant C > so that R X (X) is 
analytic in the region 

{C < ImA < Aflog|A|, |Re A| > C M » 1} (1-8) 
for every M 3> 1. Furthermore, there are infinitely many resonances in {0 < ImA < C}. 

In the case ^4 = 0, Stefanov-Vodev [12] showed that there is a free of resonances region of 
the form 

{CjvIAI - * < ImA < Mlog|A|, |Re A| > C M > 1} 

for every M, N 3> 1, while in {0 < Im A < C^vj Aj"^ } there are infinitely many resonances (called 
Rayleigh resonances) due to the Rayleigh surface waves. Later on Sjostrand-Vodev [10] proved 
that the counting function of these resonances is 

T 2 c^ 2 Vol(i> 2 + 0(r), r > 1, (1.9) 

where cr > is the speed of the Rayleigh waves and 

t 2 = (2vr)~ 2 Vol({x e R 3 : |x| = 1}). 

We expect that the counting function of the resonances in {0 < ImA < C} in the general case 
(i.e. for non-trivial A) satisfies (1.9) as well with possibly an worse bound for the remainder 
term. 

On the other hand, extending a previous result by Burq [4] to the elastic system, Bellassoued 
[3] obtained a free of resonances region of the form {0 < ImA < e _c, ' A '}, C > 0, so the Rayleigh 
resonances are concentrated in a region of the form { e - c|A| < ImA < C N \X\~ N }. Moreover, if 
the boundary T is analytic, Vodev p3] improved this region to {e~ c l A l < ImA < e- c '\ x \}. The 
presence of a non-trivial dissipative term A, however, changes the distribution of the resonances 
considerably. 

As a consequence of (1.7) we get a decay rate of the local energy of the solutions to (1.3). 

Corollary 1.2 Under the assumptions (1.1), (1-4) an d (1-5), for every a > 1, m > 0, there 
exists a constant C = C(a,m) > so that we have (for t^>l) 

\\V x u(t, -)IU 2 (n a) + \\d t u(t,-)\\ L 2 ina) < C (rtogi)™ (\\Vh\\Hm {a) + ||/ 2 ||H-(fi)) , (1-10) 

where f2 a := O n {|x| < a} and supp/j C £l a , j = 1> 2. 

The fact that (1.7) implies (1-10) was proved in [9] in the case of a unitary group. In our 
case this can be done following the approach developed in [8] (and also in [4]). Note that in the 
case ^4 = 0, Bellassoued [3] proved (1.10) with i _1 \ogt replaced by (logt) -1 . 

It turns out that if the dissipation on the boundary is stronger, we have a uniform exponential 
decay of the local energy. Indeed, consider the following mixed boundary value problem 

(d? - A e )u = in (0, +oo) x Q, 
< Bu + Ad t u = on (0, +oo) x T, (1.11) 
u(0) = fi, d t u(0) = f 2 , 

V 
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where A is as above. The outgoing resolvent, R(X), corresponding to the problem (1.11) is 
defined via the equation 

' (A e + X 2 )R(X)f = f in Q, 

(B — iXA) R(X)f = on T, (1.12) 
R(X)f — A — outgoing. 
We have the following 

Theorem 1.3 Assume (1.1) and (1.4) fulfilled with T$ = V. Then, R X (X) extends analytically 
to { | Im A | < Ci, |ReA| > C2 > 0} and satisfies there the estimate 



R X (X) 



<c'\x\- 



(1.13) 



As a consequence of (1.13) we get an exponential decay of the local energy of the solutions 
to (1.11). 

Corollary 1.4 Under the assumptions of Theorem 1.3, for every a> 1, there exist constants 
C = C(a) > 0, a > 0, so that we have (for t S> 1) 

- at l^r '-ii \ (1.14) 



\\V x u(tr)\\mn a ) + \\d t u(t,-)\\ L2{Qa) < Ce~ at (||VA|| i2(n) + \\h\\ L 2 {n) 
provided supp/j C £l a , J ' = 1> 2. 

The fact that (1.13) implies (1.14) is more or less well known in the case of unitary groups 
(e.g. see [15 ). In the case of semi-groups the proof goes in the same way (see [7]). 

It is worth noticing that an interior dissipation of the elastic wave equation with Neumann 
boundary conditions does not improve the decay of the local energy. Indeed, consider the 
following mixed boundary value problem 



(df - A e + A{x)d t )u = in (0, +00) x Q, 
Bu = on (0,+oo) x T, 
u(0) = h, d t u(0) = h, 



(1.15) 



where A G Cq°(0) is a 3 x 3 matrix- valued function satisfying the properties A = A*, A > 0. 
Then, the quasi-modes constructed in [12], [13], which are due to the existence of the Rayleigh 
waves and hence supported in an arbitrary small neighbourhood of the boundary, are also quasi- 
modes for the problem with non-trivial A. Therefore, in the same way as in these papers one 
can show that there exists an infinite sequence {Xj} with < ImAj < Cn\Xj\~ n , V7V S> 1, so 
that the following problem has a non-trivial solution: 



(A e - iXjA(x) + X])vj 
Bvj = on r, 
Xj — outgoing. 



in O, 



(1.16) 



Note finally that the situation is completely different for the usual scalar-valued wave equation 
with dissipative boundary conditions like those above. Indeed, in this case if the obstacle is 
non-trapping, the corresponding cut-off resolvent extends analytically through the real axis to 
a strip and as a consequence we have an exponential decay of the local energy without extra 
assumptions (e.g. see PQ). In other words, the behaviour of the cut-off resolvent and the local 
energy is the same as in the case of the self-adjoint problem with Neumann boundary conditions. 
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2 Proof of Theorem 1.1 

It sufices to prove (1.7) for real A S> 1, only. Let v € L^ omp (0) and let it be the solution to the 
equation 

(A e + X 2 )u = v in tt, 

(B-iA)u = on T, (2.1) 
u — A — outgoing. 

Clearly, (1.7) is equivalent to the estimate 

IMlL 2 (Q a ) < C a \\v\\ L 2(ty, A > A , (2.2) 

for every a> 1 with constants C a , Ao > indpendent of A. To prove (2.2) we need a priori 
estimates of the solutions to the equation 

(A e + A 2 )u = v in Q, 

u\ r = f,\- 1 Bu\ r = g, (2.3) 
u — A — outgoing. 

where v € L 2 omp (Q). We have the following 

Proposition 2.1 There exist constants C,Xq > so that for A > Aq we have 



u 



HHa a ) + IMLe^d < CA L \\v\\ L 2 {n) + C\\f\\ H i (r) . (2.4) 



Hereafter the Sobolev spaces H 1 are equipped with the semi-classical norm (with a small param- 
eter X^ 1 )- 

Proof. In the case of the Euclidean Laplacian A the a priori estimate (2.4) is proved in [5] 
(see Theorem 3.1). In our case the proof goes in the same way, but we will sketch it for the sake 
of completeness. Observe first that the solution to (2.3) is of the form 

u = G(\)v + K(X)f, 

where G(X)v solves the problem 

(A e + X 2 )G(X)v = v in U, 

G(X)v\ r = 0, (2.5) 
G(X)v — A — outgoing, 



while K(X)f solves the problem 



(A e + X 2 )K{X)f = in Q, 

K(X)f\ r = /, (2.6) 
K(X)f — A — outgoing. 
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Since the strictly convex obstacles are non-trapping for the Dirichlet problem of the elastic wave 
equation (see [H]), we have the estimate 

\\G(X)v\\ m{na) ^CaX-^lvWvw, A>A . (2.7) 

Thus, to prove (2.4) we need the estimate 

\\K(X)f\\ Hl{na) < C«||/|| H i (r) , A>A . (2.8) 

This in turn follows from the fact that, since the obstacle is strictly convex, one can construct 
a parametrix of K(X) near the boundary, which satisfies (2.8). More precisely, there exist a 
neighbourhood O' C Q of T and operators 

/C(A) = 0(1) : H 1 ^) -> H\Q'), K(X) = 0(A-°°) : H\T) -> H l {Sl'), (2.9) 

solving the equation 

' (A e + A 2 )/C(A)/ = ^(A)/ in ft', 
K(A)/| r = /. 

Note that such operators are constructed in [12] (Section 2). Let ip € C°°(f2), supp^ C Si', 
-0 = 1 near I\ We have 

(A e + A 2 )V>/C(A)/ = [A e ,#C(A)/ + ^(A)/ in O, 
^ ^/C(A)/| r = /, 
which leads to 

(A e + A 2 )(if(A)/-«A)/) = -[A e ,^]/C(A)/-V^(A)/ in ft, 
^ (#(A)/ - ^lC(X)f) | r = 0. 

Hence 

K (A)/ = V/C(A)/ - G(A) ([A e , V]/C(A)/ + ^(A)/) . (2.11) 

Thus, (2.8) follows from (2.11), (2.7) and (2.9). To complete the proof of (2.4) we need to show 
that 

lbllL2(r) < CA _1 ||t;|| L 2 {n) + ||u|| H i(fi a ) + C\\f\\ H i(r). (2.12) 

To this end, we write the operator A e in normal coordinates y = (yi,y r ) S R + x V in a 
neighbourhood of the boundary. We have 

dxj = Vj{y')d yi + (3j(y) ■ V y >, j = 1, 2, 3, 

where v(y') = (vi(y f ), U2(y'), 1^3 (y')) is the unit normal at y' € T. Hence 

A e = A{y')d 2 yi + Q(y, fy) + Q^y, d y ), 

where Q and Qi are second and first order differential operators, respectively, while A(y') is a 
symmetric matrix-valued function defined by 

3 

(A(y')u) k = Hou k + (X + fi )v k ^2 vjUj , k = 1,2,3. 

3=1 



6 



It is easy to check that 

det A(y') = fi 2 {\ + 2^ ) > 0. 

Set 

E(yi) = ((Q + X 2 )^u)( yil •), (ifm)( yi ,-)) L2 + (Ad yi ftMfoi, •), d yi ^u)(y u ■)) L , , 
ip being the function above. We have 

^M = ([^ 1 ,Q](^)(y 1 ,-),(^)(yi,-)) L2 

+2Re ((Q + X 2 )^u)(yi, •), d yi (ifm)(y u •)) ^ + 2Re •), ^ frMfoi. ")) L2 

= ([^i/i , Q] (i/i > •)> OH (yi > -)) L 2 - 2Re (Qi (^u) (yi , •), d yi (ijm) {yi ,-))& 
+2Re ((A e + X 2 )(iPu)( yi , •), ^foMfoi, ■)) 



' L 2 

Hence 

^(0) = - jH < 0(A 2 )||^*ll^( n) + 0(l)ll(A e + A 2 )(^)||l 2(n) . 

On the other hand 

||0 yi ftM(O, Olll^) < ^(0) + O(A 2 )||(Vu)(0, OHl^r), 
with a constant C > 0. Combining these estimates we get 

A- 1 11^(0, 011^^) < O(l)||«(0, Oliver) + 0(l)||«||Hi ( n ) + (A _1 )||(A e + A 2 H| i2(n) , 
which clearly implies (2.12). □ 

Set w = G(X)v, where v is as in (2.1). If u is the solution to (2.1), then the function u — w 
solves the equation 

(A e + X 2 )(u-w) = in n, 

(B — iA) (u — w) = —Bw on T, (2.13) 
(u — w) — X — outgoing. 

Set f = u\ r = (u- w)\ r , g = -X^Bw\ r . By (2.4), 

\\g\\mr) < C7A- 1 ||v|U a(n) , A>A . (2.14) 

Furthermore, we have 

X- 1 B(u-w)\ r = N(X)f, (2.15) 

where N(X) : H 1 ^) — > L 2 (r) is the outgoing Neumann operator. Thus, we get that the function 
/ satisfies the equation 

(N{\) - iX^A) f = g (2.16) 

with g satisfying (2.14). It is easy to see that (2.2) follows from combining (2.4), (2.14) and the 
following 



Theorem 2.2 Under the assumptions (1.1), (1-4) and (1-5), there exist constants C, Xq > so 
that the solution to (2.16) satisfies the estimate 



\\f\\m(T) < C^IMIz, 2 (r)> ^ > -V 
Proof. Since the outgoing Neumann operator satisfies 

-Im(iV(A)/,/) L 2 (r) >0, 

we obtain 



(Af,f) L 2 (r) < -Xlm{g,f) L2{T) < P~ 2 X 2 \\g 



+ P ll/llL 2 (r)' 



for every > 0. By (1.4) and (2.19), 

ll/llL 2 (r )< C(3- 1 m\L2(r) + PWfWmr)- 



(2.17) 

(2.18) 
(2.19) 

(2.20) 



Now, using (1.5) together with the properties of the outgoing Neumann operator, we will prove 
the estimate 

||/|| H i (r) <CX\\g\\ L 2 ir) + C\\f\\ L 2 (Fo) . (2.21) 

Clearly, (2.17) follows from (2.21) and (2.20) provided (3 is taken small enough. 

To prove (2.21) we will make use of the properties of the parametrix, M(X), of N(X) con- 
structed in Section 3 of [12] . First of all, we have 



||iY(A)/-A/-(A)/|| L2(r) <0(A- 



L2(D- 



(2.22) 



Moreover, M(X) is a A — ^fDO with a characteristic variety S = {( £ T*T 



-R 



} 



belonging to the elliptic region of the corresponding boundary value problem. In the region 
{C £ T*T : HCll > c^ 1 } the operator J\f{X) is an elliptic A — ^DO of class L ' (T) (hereafter we 
use the same notations as in the appendix of |12j). while in the region {( G T*T : \\(\\ < c^ 1 } it 
is hypoelliptic. Clearly, so is the operator 7V(A) — iX^A. Therefore, if x € Co°(T*r), % = 1 on 



{C G T*T : 



< e}, X = on {C G T*T : 



R 



|0 PA (1 - x)f\\m { r) < 0(A 2/3 ) J (M(X) - iX-'A) f 



> 2e}, < e < 1, we have 
i2(r) +0(A- 0O )ll/ll^(r) 



<0(A 2 / 3 )||, 



lffllL 2 (r) +0(X~ " )\\j\\m(r)- 

On the other hand, near £ the operator A^(A) is a A — ^DO of class Lp^r), whose principal 
symbol is a symmetric 3x3 matrix-valued function with eigenvalues a\(() = 2i(C)(cr||C|| — 1), 
5i(£) > 0, 02(C) > 0) °3(C) > near S. It is shown in [llj (Theorem 3.1) that there exist elliptic 
A - ^DOs, [7(A) and V(X), of class o( r ), so that we have 



(2.23) 



U(X)*M(X)Op x (x)U(X) 



-A- 2 c|A r 



A- 1 ^ - 1 








Op A (xi) mod L° ' Q 2 (T), 



(2.24) 

where — Ap denotes the positive Laplace-Beltrami operator on T, oq is a classical (independent 
of A) zero order ^DO on T with a real-valued principal symbol, and x,Xi € Cg°(T*r), xi = 1 



on {C G T*r : 



R 



< 3e}, xi = on {C G T*T 



R 



> 4e}, x = 1 on suppxi- In 
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fact in [TT] a better diagonalization of A/*(A) near S is carried out, but for our purposes (2.24) 
will suffice. Now the function / = U(X)~ 1 Opx(x)f satisfies 



-A~ 2 c|A r - A _1 a - 1 
V{X) 



f-i\~ l Af = g, (2.25) 



where A = U*AU is a A — ^DO of class Lq'q(F) with a principal symbol satisfying <r p (A) > 0, 
tjp(A)* = cr p (A), and g satisfies 

llffll^(r) < Ibll^(r) + 0(A-°°)||/|| H i (r) . (2.26) 

Writing / = (fufafa), 9 = (91,92,93), we reduce (2.25) to 

(-A- 2 c|A r - 1 - iX-%) fi=9i + ^ 1 (ba/a + 63/3) , (2-27) 

I:)- 1 (::;:)• 

where fy, Cj are scalar- valued A — ^>DOs of class L ' (r), the principal symbol of b\ satisfying 
Rec7p(6i) > 0, while V(X) is a 2 x 2 matrix-valued elliptic A — ^DO of class Lqq(F). Thus, the 
inverse V(X)~ 1 is again a 2 x 2 matrix-valued elliptic A — ^DO of class L Q ' (T), so we can solve 
the equation (2.28). In particular, we obtain 

ll/2||L 2 (r) + Whh^r) < Il52|k 2 (r) + ll53||L 2 (r) + 0(X~ 1 )\\fi\\ L 2 ( py ( 2 -29) 

Furthermore, we conclude that the function f\ solves an equation of the form 

-A~ 2 4A r - 1 - iX~ l b\ /1 = h, (2.30) 

where b is a scalar- valued A — ^>DO of class Lq 'g(r) with a principal symbol satisfying Re <r p (b) > 
0, and h satisfies 

\\h\\mn < hh^D + o(x-™)\\f\\ m{r) . (2.31) 

We are going to show that the assumption (1.5) leads to the estimate 

||/i||L 2 (r) < CX\\h\\ L 2 (r) + C||/ 1 || i2(n) . (2.32) 

Before doing so, observe that (2.32) implies (2.21). Indeed, since 

ll/i||L 2 (r' )<C||/lli 2 (r )+O(A- oo )||/|| i2(r) , 

we deduce from (2.29), (2.31) and (2.32) that 

l|Op A (x)/||^ ( r) < CX\\g\\ LHr) + C||/|| x2(ro) + 0(X-°°)\\f\\ m{T) . (2.33) 

Thus, (2.21) follows from (2.23) and (2.33). 

The fact that (1.5) implies (2.32) can be derived from the more general results of [2], but 
we will give here a simpler proof following [16] where this is carried out in the case b = (see 
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Theorem 2.3 of [IS]). Denote by r (x,£), G T*T, the principal symbol of the operator 

— cfjAr, so that we have £ = {(x,£) G T*T : ro(x,£) = 1}. Recall that the bicharacteristic flow 
$(t) : T*T -> T*T, t G R, associated to the Hamiltonian r (x,£) is defined by $(t)(x°,^°) := 
(x(i), £(i)), where the pair (x(t),£(t)) solves the Hamilton equation 



dx(t) 9ro(x,^) d£(t) dro(x,£) 



Of 



Of, 



dx 



, x(o) = x o ,m = a°- 



Fix a point £ c 



n0 ,£ ) G S and choose a real-valued function p(x,£) G C^°(T*r), < p < 1, 



such that p = 1 in a neighbourhood of C° and p = outside a biger neighbourhood. Given a 
t £ R, define the function P t{x,£) G Cq°(T*T) by Pt{x,Cj = p(<&(— f)(:r, £)). By a microlocal 
partition of the unity in a neighbourhood of E, it is easy to see that (2.32) follows from (1.5) 
and the following 

Lemma 2.3 For every T > there exist positive constants C = C(T) and Xq = \q(T) so that 
the solutions to (2.30) satisfy the estimate 



p(x,V x )fi 



< 



p t (x,V x )fi 



L2(r) 



+ 2T\\\h\\ L2{T) + C\- 1 h 



l"(t) 



(2.34) 



for < t < T , X > Xq. Hereafter we denote T> x := X l D x . 
Proof. Set P = -X~ 2 c 2 R A r - 1 - iA -1 &. Since 

dtPt + {ro,Pt} = 0, 

the operator 

Qt := XdtPt{x,V x ) +iX 2 [P,p t (x,V x )] 

is a zero order A — ^DO, and hence uniformly bounded on L 2 (T). Moreover, the fact that the 
principal symbol of the operator b satisfies Kea p (b) > implies 

-Re(6/,/) L2(r) < 0(A- 1 )||/||| 2(r) , V/ G L 2 (Y). 

Therefore, using the identity 

\j t \\pt{x,v x) r 2 



l ., =Re {dtPt{x,V x )fi,p t (x,V x )fi / l ^ 



ARe ([P, Pt (x,V x )]f 1 ,p t (x,V x )f 1 ) L r) + A- 1 Re {QtfuPtix, V x )f 1/ ^ | 



ARe (pt(x,V x )Ph,pt(x,V x )f; 



L2(r) 

Re (bpt{x,T> x )fi,pt(x,V x )fi 



+ A~ 1 Re (Qtfi,Pt{x,V x )h 



we obtain 



By (2.35), 



dt 



Pt{x,V x )fi 



L2(r) 



< 2A 



Pfi 



L"(T) 



+ 0(X- 1 ) h 



(2.35) 



L»(r) 



d_ 
~d7 



Pr{x,V x )fi 



dr 
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< 



p t (x,V x )fi 



L 2 (r) 



+ 2tX 



Pfi 



L2(r) 



+ 0(X~ l ) h 



L2(r) 



□ 



Clearly, the free of resonances region follows from the following 



Proposition 2.4 Under the assumption (1-1), for X belonging to the region (1-8) with a suitably 
chosen constant C > 0, the solution to (2.16) satisfies the estimate 



a 



(2.36) 



L2 ^ S hnT miL2 ^ 

Proof. Without loss of generality we may suppose Ai := Re A 3> 1. It is shown in [12] that, 
for A belonging to the region Am = {0 < ImA < Mlog |A|, Re A > Cm ^> 1}, the Neumann 
operator has a parametrix Af(X) which is a Ai — *$>DO with a characteristic variety S, depending 
on a parameter AT^Im A <C 1. In particular, (2.22) still holds with 0(X~°°) replaced by 0(A^°°). 
In what follows we will keep the same notations as in the proof of Theorem 2.2 above. In fact, 
much of the analysis still works with Op A replaced by Op Al . For example, we have the following 
analogue of (2.23) 



|0 PAl (l - x)/ll^cr) ^ 0(|A| 2 / 3 )|| 5 || L2(r) + 0(Ar 



tfi(r), 



A € A 



M- 



(2.37) 



We still have (2.29) with 0(A _1 ) replaced by C^Af 1 ) as well as (2.30) with h satisfying (2.31) 
with 0(A-°°) replaced by 0(X^°°). Thus, we get 



ImAiMli^ = -Im + ( Afc/ : . , 

< l^| 2 ||^llL2(r)ll/ilk 2 (r) + C|^lll^ll|2(r)> 



and hence 



ImAHMU^D ^ (2ImA-C)||/i|| L 2 (r) < 0(|A|)||/i|| L2(r) , 
provided ImA > C, A 6 Am- Combining (2.38) with (2.29) and (2.31), we obtain 

C"| 



(2.38) 



l|Op Al (x)/|| L2( r) < j^\\g\\ L Hr) + 0(X^)\\f\\ HHr) , 



for A belonging to the region (1.8). Now (2.36) follows from (2.37) and (2.39). 

To prove the existence of infinitely many resonances (i.e. poles of (AA^(A) 
{0 < Im A < C} we will proceed as in [12]. Without loss of generality we may suppose Re A > 
By (2.36), we have 



(2.39) 
iA)' 1 ) in 



(AiV(A) - iA)- 1 



L 2 (r)^L 2 (r) 



<C(log|A|)-\ AG/i, 



(2.40) 



where := {A € C : ±ImA = log Re A, Re A > C'} with some constant C S> 1. If we 
suppose that (XN(X) — iA)' 1 is analytic in {A € C : < Im A < C, Re A > C'}, so it is in 
{A G C : |ImA| < log Re A, Re A > C'}. Then, by (2.40) together with the Fragmen-Lindelof 
principle we get 



(AiV(A) - iA)' 1 



L 2 (r)^L 2 (r) 



<C(logA)-\ AGR, X>C. 



(2.41) 
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On the other hand, it is shown in [T2] that there exist quasi-modes (fj,kj) G L 2 (T) x R such 
that ||/j||/,2 = 1, kj — » +00 and 

\\kjN(kj)fj \\ L 2 < Const. 

Hence, 

WikjNikj) -iA)fj\\ L2 < Const, 

which combined with (2.41) lead to 

l = \\f j \\ L 2<Const{\ogk j )-\ 

which is impossible if we take kj large enough. Therefore, the operator-valued function (AiV(A) — 
iA)~ l cannot be analytic in {A G C : < Im A < C, Re A > C'}. □ 

3 Proof of Theorem 1.3 

Again, it suffices to prove (1.13) for real A S> 1, only. Let v G L^ omp (0) and let u be the solution 
to the equation 

(A e + X 2 )u = v in 

{B — iXA) u = on T, (3.1) 
u — A — outgoing. 
Clearly, (1.13) is equivalent to the estimate 

IM|x2(n a ) < CaA" 1 11^11^2(0) , A > A . (3.2) 
The function f = u\r solves the equation 

(N(X) -iA)f = g (3.3) 

with g satisfying 

IMIz^r) < CX^WvW^^, X > A . (3.4) 
Thus, in view of (2.4), to prove (3.2) it suffices to show that 

WfWHHF) <C\\g\\ L 2 {r) , A>A . (3.5) 
Using (2.18) and (1.4) with T = T, we get 

C||/||| 2{r) < (Af,f) LHr) < -lm(gJ) L 2 {r) <r 2 \\g\\h (r) +(3 2 \\f\\h ir) , (3.6) 
for every f3 > 0. Taking (3 small enough, we deduce from (3.6), 

WfWmv) < C\\g\\ L 2 {r) . (3.7) 

Let 77 G C^(T*T), n = 1 on {C G T*r : ||C|| < c^ 1 }, ?] = on {( S T*T : ||C|| > c^ 1 + 2}. Since 
the parametrix M(X) on supp (1 — rj) is an elliptic A — ^DO of class L 'o(r), we have 

[|Op A (l - v)f\\m { r) < C \\N(X)f\\ L 2 {T) + 0{X-°°)\\f\\ m{T) 

< C\\g\\ L 2 [r) + C||/|| L2( r) + 0(X-°°)\\f\\ HHr) . (3.8) 

On the other hand, 

l|Op A (r7)/||^ (r) <C||/|| x2(r) . (3.9) 

Now (3.5) follows from combining (3.8) and (3.9) with (3.7). □ 
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